Search for: All records

Assigning gene function from genome sequences is a rate-limiting step in molecular biology research. A protein's position within an interaction network can potentially provide insights into its molecular mechanisms. Phylogenetic analysis of evolutionary rate covariation (ERC) in protein sequence has been shown to be effective for large-scale prediction of functional relationships and interactions. However, gene duplication, gene loss, and other sources of phylogenetic incongruence are barriers for analyzing ERC on a genome-wide basis. Here, we developed ERCnet, a bioinformatic program designed to overcome these challenges, facilitating efficient all-versus-all ERC analyses for large protein sequence datasets. We simulated proteome datasets and found that ERCnet achieves combined false positive and negative error rates well below 10% and that our novel “branch-by-branch” length measurements outperforms “root-to-tip” approaches in most cases, offering a valuable new strategy for performing ERC. We also compiled a sample set of 35 angiosperm genomes to test the performance of ERCnet on empirical data, including its sensitivity to user-defined analysis parameters such as input dataset size and branch-length measurement strategy. We investigated the overlap between ERCnet runs with different species samples to understand how species number and composition affect predicted interactions and to identify the protein sets that consistently exhibit ERC across angiosperms. Our systematic exploration of the performance of ERCnet provides a roadmap for design of future ERC analyses to predict functional interactions in a wide array of genomic datasets. ERCnet code is freely available at https://github.com/EvanForsythe/ERCnet. 

We hypothesize that bacteria isolated from free-ranging animals could potentially be useful for practical applications. To meet this objective a Gram-positive bacterium was isolated from the gastrointestinal (GI) tract of a Gray Wolf (Canis lupus) using Brucella broth with hemin and vitamin K (BBHK). By small ribosomal RNA (16S) gene sequencing the bacterium was initially identified as a novel Carnobacterium maltaromaticum strain. The bacterium could be propagated both anaerobically and aerobically and was both catalase/oxidase negative and negative by the starch hydrolysis as well as negative using lipase assays. The reference whole genome sequence (WGS) was obtained using both Illumina and Nanopore sequencing. The genome assembly was 3,512,202 bp in length, encoding core bacterial genes with a GC% content of 34.48. No lysogenic bacteriophage genes were detected, although the genome harbors genes for the expression of bacteriocin and other secondary metabolites with potential antimicrobial properties. Multilocus sequence typing (MLST), WGS phylogenetics, average nucleotide identity (ANI), and single nucleotide polymorphism (SNP) analyses of the isolate’s genome indicate this bacterium is a newly identified Carnobacterium maltaromaticum sequence type (ST). Members of the Carnobacteria have anti-listeria activities, highlighting their potential functional properties. Consequently, the isolate could be a potential probiotic for canids and this is the first report on an axenic C. maltaromaticum culture from the genus Canis. 

Inflammatory bowel disease (IBD) is increasing among mammals around the world, and domestic dogs are no exception. There is no approved cure for canine IBD with limited treatment options. Novel probiotic bacteria discovery from free-ranging animals for the treatment of IBD in domestic pets can likely yield promising probiotic candidates. Consequently, the overall aim was to isolate bacteria from free-ranging animals that could potentially be utilized as novel probiotics. Two bacteria identified as unique Paenibacillus spp. strains by small ribosomal RNA (16S) gene sequencing were isolated from the gastrointestinal tract of a North American Gray Wolf (Canis lupus). The bacteria were typed as Gram-variable, and both were catalase/oxidase positive as well as sensitive to commonly used antibiotics. The bacteria digested complex carbohydrates and lipids by standard assays. The isolated bacteria also inhibited the growth of Staphylococcus aureus and Micrococcus luteus. The whole genome sequence (WGS) length of bacterial isolate ClWae17B was 6,939,193 bp, while ClWae19 was 7,032,512 bp, both similar in size to other Paenibacillus spp. The genomes of both bacteria encoded enzymes involved with the metabolism of complex starches and lipids, such as lyases and pectinases, along with encoding antimicrobials such as lanthipeptides, lasso peptides, and cyclic-lactone-autoinducers. No pernicious virulence genes were identified in the WGS of either bacterial isolate. Phylogenetically, the most closely related bacteria based on 16S gene sequences and WGS were P. taichungensis for ClWae17B and P. amylolyticus for ClWae19. WGS analyses and phenotypic assays supported the hypothesis that the isolates described constitute two novel candidate probiotic bacteria for potential use in dogs. 

The interaction between the nuclear and chloroplast genomes in plants is crucial for preserving essential cellular functions in the face of varying rates of mutation, levels of selection, and modes of transmission. Despite this, identifying nuclear genes that coevolve with chloroplast genomes at a genome-wide level has remained a challenge. In this study, we conducted an evolutionary rate covariation analysis to identify candidate nuclear genes coevolving with chloroplast genomes in Juglandaceae. Our analysis was based on 4,894 orthologous nuclear genes and 76 genes across seven chloroplast partitions in nine Juglandaceae species. Our results indicated that 1,369 (27.97%) of the nuclear genes demonstrated signatures of coevolution, with the Ycf1/2 partition yielding the largest number of hits (765) and the ClpP1 partition yielding the fewest (13). These hits were found to be significantly enriched in biological processes related to leaf development, photoperiodism, and response to abiotic stress. Among the seven partitions, AccD, ClpP1, MatK, and RNA polymerase partitions and their respective hits exhibited a narrow range, characterized by dN/dS values below 1. In contrast, the Ribosomal, Photosynthesis, Ycf1/2 partitions and their corresponding hits, displayed a broader range of dN/dS values, with certain values exceeding 1. Our findings highlight the differences in the number of candidate nuclear genes coevolving with the seven chloroplast partitions in Juglandaceae species and the correlation between the evolution rates of these genes and their corresponding chloroplast partitions. 

Abstract There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes—even single-nucleotide substitutions—create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein–protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution. 

The discovery of novel probiotic bacteria from free-ranging animals for the treatment of inflammatory bowel disease in domestic pets is a unique approach. The chloroform extraction of gastrointestinal (GI) tract material was used to inactivate vegetative cells and select for spore-forming bacteria. A bacterium identified as a novel Paenibacillus sp. strain via small ribosomal RNA (16S) gene sequencing was isolated from the GI tract of a gray wolf (Canis lupus). The bacterium was typed as Gram-variable, both catalase/oxidase-positive and positive via starch hydrolysis and lipase assays. The bacterium inhibited the growth of Staphylococcus aureus, Escherichia coli and Micrococcus luteus. The draft whole genome sequence (WGS) assembly was 7,034,206 bp in length, encoding 6543 genes, and is similar in size and coding capacity to other closely related Paenibacillus spp. The isolate’s genome encodes several germination and sporulation gene products along with antimicrobials such as a bacteriocin system and chitinase. Enzyme genes such as alpha amylase, cellulase, lipases and pectin lyase are also present in the genome. An incomplete lysogenic bacteriophage genome was also present in the isolate’s genome. Phenotypic characteristics combined with a WGS genotype analysis indicate that this bacterium, designated Paenibacillus sp. ClWae2A, could be a potential candidate probiotic for domestic dogs. 

Mitochondrial and plastid functions depend on coordinated expression of proteins encoded by genomic compartments that have radical differences in copy number of organellar and nuclear genomes. In polyploids, doubling of the nuclear genome may add challenges to maintaining balanced expression of proteins involved in cytonuclear interactions. Here, we use ribo-depleted RNA sequencing (RNA-seq) to analyze transcript abundance for nuclear and organellar genomes in leaf tissue from four different polyploid angiosperms and their close diploid relatives. We find that even though plastid genomes contain <1% of the number of genes in the nuclear genome, they generate the majority (69.9 to 82.3%) of messenger RNA (mRNA) transcripts in the cell. Mitochondrial genes are responsible for a much smaller percentage (1.3 to 3.7%) of the leaf mRNA pool but still produce much higher transcript abundances per gene compared to nuclear genome. Nuclear genes encoding proteins that functionally interact with mitochondrial or plastid gene products exhibit mRNA expression levels that are consistently more than 10-fold lower than their organellar counterparts, indicating an extreme cytonuclear imbalance at the RNA level despite the predominance of equimolar interactions at the protein level. Nevertheless, interacting nuclear and organellar genes show strongly correlated transcript abundances across functional categories, suggesting that the observed mRNA stoichiometric imbalance does not preclude coordination of cytonuclear expression. Finally, we show that nuclear genome doubling does not alter the cytonuclear expression ratios observed in diploid relatives in consistent or systematic ways, indicating that successful polyploid plants are able to compensate for cytonuclear perturbations associated with nuclear genome doubling. 

Abstract Cytonuclear coevolution is a common feature among plants, which coordinates gene expression and protein products between the nucleus and organelles. Consequently, lineage-specific differences may result in incompatibilities between the nucleus and cytoplasm in hybrid taxa. Allopolyploidy is also a common phenomenon in plant evolution. The hybrid nature of allopolyploids may result in cytonuclear incompatibilities, but the massive nuclear redundancy created during polyploidy affords additional avenues for resolving cytonuclear conflict (i.e. cytonuclear accommodation). Here we evaluate expression changes in organelle-targeted nuclear genes for 6 allopolyploid lineages that represent 4 genera (i.e. Arabidopsis, Arachis, Chenopodium, and Gossypium) and encompass a range in polyploid ages. Because incompatibilities between the nucleus and cytoplasm could potentially result in biases toward the maternal homoeolog and/or maternal expression level, we evaluate patterns of homoeolog usage, expression bias, and expression-level dominance in cytonuclear genes relative to the background of noncytonuclear expression changes and to the diploid parents. Although we find subsets of cytonuclear genes in most lineages that match our expectations of maternal preference, these observations are not consistent among either allopolyploids or categories of organelle-targeted genes. Our results indicate that cytonuclear expression evolution may be subtle and variable among genera and genes, likely reflecting a diversity of mechanisms to resolve nuclear-cytoplasmic incompatibilities in allopolyploid species. 

Abstract Nuclear and plastid (chloroplast) genomes experience different mutation rates, levels of selection, and transmission modes, yet key cellular functions depend on their coordinated interactions. Functionally related proteins often show correlated changes in rates of sequence evolution across a phylogeny (evolutionary rate covariation or ERC), offering a means to detect previously unidentified suites of coevolving and cofunctional genes. We performed phylogenomic analyses across angiosperm diversity, scanning the nuclear genome for genes that exhibit ERC with plastid genes. As expected, the strongest hits were highly enriched for genes encoding plastid-targeted proteins, providing evidence that cytonuclear interactions affect rates of molecular evolution at genome-wide scales. Many identified nuclear genes functioned in post-transcriptional regulation and the maintenance of protein homeostasis (proteostasis), including protein translation (in both the plastid and cytosol), import, quality control and turnover. We also identified nuclear genes that exhibit strong signatures of coevolution with the plastid genome, but their encoded proteins lack organellar-targeting annotations, making them candidates for having previously undescribed roles in plastids. In sum, our genome-wide analyses reveal that plastid-nuclear coevolution extends beyond the intimate molecular interactions within chloroplast enzyme complexes and may be driven by frequent rewiring of the machinery responsible for maintenance of plastid proteostasis in angiosperms.